Symmetry of the spin Hamiltonian for herbertsmithite: a spin-| kagome lattice 
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We present magnetization measurements on oriented powder of ZnCu3(OH)6Cl2 along and per- 
pendicular to the orienting field. We find a dramatic difference in the magnetization between the 
two directions. It is biggest at low measurement fields H or high temperatures. We show that the 
difference at high temperatures must emerge from Ising-like exchange anisotropy. This allows us to 
explain muon spin rotation data at T — > in terms of an exotic ferromagnetic ground state. 

PACS numbers: 75.50.Lk, 75.10.Nr 



The synthesis of the herbertsmithite [ZnCu 3 (OH) 6 Cl2] 
[l| has led to a renewed interest in the frustrated spin- 
1/2 Heisenberg model on the kagome lattice. This sys- 
tem has a highly degenerate ground state Q and any 
small perturbation to the Hamiltonian can severely af- 
fect the ground state manifold. The perturbations can 
be: exchange anisotropy 0, bond anisotropy 0, 
transverse field [H, 0] , Dzyaloshinksy-Moriya Interaction 
(DMI) mo, or longer range interactions [ll| . There- 
fore, numerous theoretical directions have been taken 
to predict the low-temperature behavior of the kagome 
system, and some of them were particularly applied to 
magnetization and other measurements of the herbert- 
smithite 03 03- 

This mithite is exciting since Cu ions create a spin- 
1/2 magnetic kagome layer separated by non magnetic 
Zn atoms from the adjacent layers. The compound was 
found to be a quantum spin liquid with no broken con- 
tinuous symmetry but gapless excitations 03 EH EE 
flfl [l7j . At high temperatures the inverse susceptibil- 
ity obeys a Curie- Weiss (CW) law, X = C/(T + 6), 
where C is the Curie constant and the CW temperature 
= 314 K. Below ~ 75 K a sharp increase in the suscep- 
tibility occurs, deviating from the ideal kagome Heisen- 
berg model pj| . This upturn was accounted for by DMI 
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FIG. 1: (Color online) X-ray diffraction of powder (black) 
and oriented powder (gray) from a surface perpendicular to 
the orienting field. 



[Tol . E|| or anisotropy in the bonds 0, 0] • It was also sug- 
gested that impurities from a Zn/Cu substitution play 
a significant role in the low-temperature susceptibility 
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2fj, |21j. However, free impurities, or interacting im- 
purities th at g enerate an additional ferromagnetic Curie- 
Weiss law [10(|l2j, have been shown not to describe this 
upturn completely. In fact, Rietveld refinement of our 
sample showed no Zn/Cu substitution within the exper- 
imental resolution [22| . The sample is made by the same 
procedure and group as the samples in Refs. [l5| and 
[2^| . Finally, different local probes such as muon 
O, [13] Cu, and CI [H[ nuclear magnetic resonance, and 
electron spin resonance suggest different behavior of 
the susceptibility below ~ 50 K. Thus, there is still no 
agreement on the interactions that control the magnetic 
properties of herbertsmithite. 

In fact, since it is only available as a powder, the sym- 
metries of its spin Hamiltonian are not known. To clar- 
ify these symmetries we present magnetization measure- 
ments on oriented powder of ZnCu3(OH)6Ci2 along (£) 
and perpendicular to (_L) the orienting field. The sym- 
metry of the interactions are probed at high tempera- 
tures where impurities are not expected to contribute to 
the susceptibility even if they do exist, and all probes 
roughly agree. 

The orientation was done by curing ZnCu3(OH)gCl2 
powder overnight with Stycast in a field of 8 T at room 
temperature. The samples were cured in a Teflon form 
producing a ball 6 mm in diameter. During the first 40 
minutes of the orientation, a shaking mechanism was ap- 
plied to the sample form. A particularly small amount 
of powder was used to avoid saturating the Stycast and 
eliminating powder residues at the bottom of the ball. 
We prepared a second "test" sample in the same manner, 
but this time without orientation. We refer to the second 
ball as the powder sample. The mass of ZnCu3(OH)gCl2 
in the ball is known only roughly and the absolute value 
of the molar magnetization is not accurate. We also pre- 
pared a ball made of Stycast only. All samples were mea- 
sured in a gelatin capsule. 

In Fig. [1] we plot the x-ray diffraction from the powder 
and oriented samples. For these measurements a separate 
surface perpendicular to the orienting field was prepared 
and used. The Bragg peak intensities are shown in the 
figure. In the oriented case the (002) and (006) peaks 
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FIG. 2: (Color online) Normalized magnetization M/H = x 
times the temperature versus temperature at external field 
of 400 G for oriented sample (diamonds) in two directions, 
powder sample (squares) in two directions as if it was ori- 
ented, and stycast sample (open circles). The inset shows the 
magnetization as a function of applied field for two different 
temperatures and directions for low fields. 



increased dramatically, while many of the other peaks did 
not. This x-ray picture shows a high degree of orientation 
such that the c direction is parallel to the field. The level 
of orientation will be discussed further below. 

DC magnetization measurements, M, were performed 
using a Cryogenic SQUID magnetometer in two config- 
urations. One configuration 'z' is when the orienting 
and the applied (SQUID) fields coincide, (i/||c). The 
other configuration '_L' is when the ball is rotated by 
90° and thus the applied field is in the kagome plane, 
(H _L c). In Fig. [2] we present MT/H of the two sam- 
ples, powder and oriented balls. In the reset of this 
paper we use x to indicate the normalized magnetiza- 
tion M/H (and not dM/dH). These measurements were 
taken at H = 400 G. The measurements are conducted 
as follows: we first measured the powder sample and 
then the oriented sample in both configurations. Fi- 
nally, we repeated the powder measurements for a sec- 
ond time, but rotated the powder ball as if it was ori- 
ented. All powder measurements collapse into a single 
curve, as expected, demonstrating the reproducibility of 
the measurement. The Stycast sample showed a very 
small diamagnetic signal which is also depicted in Fig. [5] 
The core diamagnetic susceptibility of ZnCu3(OH) 6 Cl2 
is — 16.7 X 10~ 5 cm 3 /mole [2J]. The Van-Vleck contribu- 
tion is expected to be of the same order of magnitude, but 
with a positive sign (2f| . Both are much smaller than the 
measured x a t room temperature of 1 x 10~ 3 cm /mole. 

In Fig. [2] no special energy scale is found in either one 
of the measurements. The only indication of an interac- 
tion between spins is the fact that xT f° r both directions 
and the powder decreases with decreasing T. \T of the 
powder is smaller than XzT and larger than x±T of the 
oriented sample. However, a comparison of the absolute 
value of x °f the powder and the oriented sample is not 
accurate. We did try to have an equal amount of sample 
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FIG. 3: (Color online) An Arrott plot showing M 2 versus 
H/M at various temperatures in the perpendicular direction 
(a) and parallel direction (b). 



in both balls but there is no telling how successful we 
were. A more relevant comparison is between x m the 
different directions of the oriented sample; % 2 T increases 
faster than x±T, and at room temperature Xz — 1-6x_l- 
Thus the ratio between the z and _L directions increases 
as the temperature increases. 

In the inset of Fig. [2] we show the magnetization as a 
function of H for two different temperatures and direc- 
tions. The magnetization is beginning to show signs of 
saturation, suggesting contribution from ferromagnetic 
impurities. To check this possibility we present in Fig. [3] 
an Arrott plot [26]. This plot takes advantage of the 
high field data. At a ferromagnetic transition, M 2 (T C ) 
is expected to be a linear function of H/M. We found 
no evidence for such linear behavior. In fact, M 2 (T) is 
independent of H/M near the origin as expected when 
the ferromagnetic critical temperature is lower than the 
available temperature. This indicates the lack of ferro- 
magnetic impurities in our measurements. 

In Fig.|4]we plot versus temperature for two fields, 
2000 and 100 G and for the two orientations. In the 
inset of Fig. HJb) we plot the at low temperatures 
(T < 50 K); clearly, Xz linearizes at T ~ 30 K whereas 
X ± linearizes at a much higher temperature (T ~ 100 K). 
6, and C in arbitrary units are extracted from a linear 
fit of the high-temperature (150 < T < 280 K) data to 
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X±\ — (T + ^i,z)/Ci !2 . The fits are shown by the solid 
line. 

In Fig. [5]we plot 9j_,z> and y/C±, z which is proportional 
to the g± tZ factor (if the sample was fully oriented) versus 
the applied field. 9± increases slowly with decreasing ap- 
plied field and saturates below 400 G. On the other hand, 
9 Z increases rapidly below 2 kG. The Curie constant has 
a similar behavior. The powder average of 9j_, z at low 
fields does not reconcile with 9 ~ 300 K measured in a 
powder and there must be some extrinsic contribution 
to the normalized magnetization in the partially aligned 
samples at low fields. However, we have no evidence that 
this contribution is due to impurities. 

In contrast, at high fields, H > 2 kG, 9 of the two 
directions is hardly distinguishable and on the order of 
the powder value. In addition, useful information can be 
extracted from the CW temperature only if it is obtained 
by measurements at T > 9. Therefore, we concentrate 
on the results obtained by H > 2 kG, as shown in the 
inset of Fig. [3 At 2 kG the ratio of ^fcTJCl = 1.179(6) 
and 9 z /9 ± = 1.23(1). 
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FIG. 4: (Color online) Inverse normalized magnetization 
H/M = x _1 versus temperature at H = 2 kG (a) and at 
H = 100 G (b). The solid lines are linear fits to the high 
temperature (> 150 K) data, (a) inset displays the inverse 
of the normalized magnetizations 1/xz (black squares) and 
VX-L fe ra y diamonds), and the inverse intrinsic normalized 
magnetization 1/xl (black line) and 1/x!l (g rav line) obtained 
from Eq. [T] The dashed lines demonstrate that 9z < 9 % z and 
8± > 8 l ± . In the inset of (b) we plot the low-temperature 
behavior of x _1 at 100 G. 
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FIG. 5: (Color online) The Currie-Wiess temperatures (filled 
symbols) and square root of the Currie constant (open sym- 
bols) of the oriented sample perpendicular to (black) and in 
the kagome plane (gray) . The inset show a zoom on the high 
field data. 



In order to convert the measured x presented above 
to the intrinsic normalized magnetization x 1 m differ- 
ent directions, it is important to estimate the level of 
orientation. This can be done using the x-ray data. 
The ratio of the x-ray intensity (I) from the powder 
R = l(00h)/ I(kk0) represents the signal intensity ra- 
tio between the two kinds of plane. Let's assume that 
there are N grains composed of two sets, aN that can 
orient perfectly with the field, and (1 — a)N that are 
not effected by the field at all since they are made of 
a few crystalline, for example. We further define (3 as 
the probability that a particular plane will contribute to 
the scattering intensity in a powder. After orientation 
the x-ray intensity ratio between the same planes would 
be R' = [aR + (3(1 - a)R\/ f3{\ - a). We can estimate 
from the width of the peaks which is 0.2° out of 180, thus 
f3 ~ 0.001. Using R and Rl of the (006) and (220) we 
find a = 0.25. This level of orientation is in agreement 
with Imai et a/.[23|. 

In an oriented sample we expect 

Xz ,± = (i - a) Qxt + ^xij + «xi,±- (1) 

This relation could be inverted to produce X z ±- in the 
inset of Fig. Ufa) we present both l/xz,X and 1/xl x f° r 
the normalized magnetization data taken at H = 2 kG. 
New intrinsic CW temperatures 9\ ± could be obtained 
from l/xl ± as demonstrated by the dashed lines. 9 l z ± 
represent the CW temperature as if the sample was fully 
oriented. Although a is just an estimate of the level 
of orientation, the important point is that 9\ > 9 Z and 

ei<9 ± . 

We now turn to discuss the possible origin of the 
X anisotropy in terms of superexchange anisotropy and 
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DMI. The DMI Hamiltonian is given by, 

H= JS i ■ S i + D « • (Si x Sj) (2) 



<%j> 



where Djj is a vector assigned to each bond. In the mean 
field approximation (Sj — > M/p/xs) this Hamiltonian is 
written as 



W = gu-B ^ ■ H 



where 



H e " = 



( JM + D x M) + H 



(3) 



D ={1/Z) Djj, and Z is the number of neighbors. 
Special attention must be taken for the convention of the 
ij bond direction since it sets the direction of Dy [Toj | . 
The magnetization satisfy the equation 

M = ?(i (JM + DXM)+H ) (4) 

where C = (gus) 2 S(S + l)/(3fcs) is the Curie constant. 
Up to first order in D 
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(T — 6 CW ) \ T — 6 CW 
where 8 CW — CZJ/(cj^,b) 2 and 

CZ 



A H. 



(5) 



In particular 
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Therefore, Dy does not contribute to the CW law. 

In contrast, the superexchange anisotropy Hamiltonian 
is given by 



(8) 



<«,j> 



In this case, if the sample was perfectly oriented, we 
would have 6 Zi ± = J z ,±/ks- Since our sample is not 
perfectly oriented, our high-temperature high field linear 
fits of X± z measures a lower bound on J z and an upper 
bound on J±. 

The lower bound on J z is larger than the upper bound 
on Ji . Despite the fact that measurements of Xz and x± 
are contaminated with x± an d xl respectively, as indi- 
cated by Eq. [IJ the conclusion J z > J± is unavoidable. It 
is robust even against possible core and Van-Vleck cor- 
rections. Thus herbertsmithite has an Ising-like exchange 
anisotropy. This, however, is not the end of the story. If 



J z > Jj_, we would expect Xz < X_L, in contrast to obser- 
vation. Therefore, to explain the high x in the z direction 
we must invoke an anisotropic g factor as well. 

In the classical ground state of antiferromagnets on 
the kagome lattice with exchange anisotropy, the spins 
are coplanar and two angles between spins (p on each tri- 
angle obey cos ip = —J z /{Jz + J±)- The third angle com- 
pletes the circle. This condition maintains the ground 
state macroscopic degeneracy. Nevertheless, unlike in 
the Heisenberg case, there is a critical temperature T c 
below which an exotic ferromagnetic order exists with fi- 
nite total magnetization, but no sublattice long-range or- 
der (27j . Upon cooling through T c the magnetization in- 
creases abru ptly and continuously down to T — > where 
it saturates [28|]. In zero field, domains can be formed, 
but a small applied magnetic field will stabilize the mo- 
ment. The powder average of the moment projection on 
the field direction is given by the value 
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6 



(1 + 2cos<^) 



(9) 



per spin. 

We believe that this ferromagnetic order contributes to 
the observed % at T — > by transverse field (TF) muon 
spin rotation (/LtSR) experiment [l4T |. In /iSR, impurities, 
if they exist, are expected to contribute to the muon line 
width while most of the sample contributes to the line 
shift. In what follows we examine what part of the /zSR 
shift can be explained by exchange anisotropy only. A 
complete understanding will of course require taking DMI 
interaction into account as well. 

The /iSK measurements were done at a field of H = 
2 kG and the shift K in the muon rotation frequency 
as a function of temperature was measured. This shift 
is a consequence of the sample magnetization; therefore, 
K is expected to be proportional to normalized magne- 
tization. The high temperature data are used to cal- 
ibrate the proportionality constant between K and x 
[l4| . The data are reproduced in Fig. [BJ \ increases 
sharply with decreasing temperatures between ~ 10 K 
and ~ 1 K and saturates below T ~ 200 mK at a value 
of x = 15.7(5) x 10 -3 cm 3 /mol Cu. This x mounts to an 
average moment of 0.006xts per Cu, in the direction of 
the applied 2 kG field. Solving Eq. [9] for the anisotropies 
gives J z /J±_ = 1.06. In Fig. [B]we present simulations de- 
scribed in Ref. [27|, for J z /J 1 _ = 1.04 and J z /J x = 1.08 
showing similar behaviour as the experiment. For this 
type of exchange anisotropy the expected T c / J± = 0.03 
as shown in the inset of Fig. [S]also taken from Ref. [27| . 
For J± ~ 200 K we expect T c = 6 K. This temperature 
is at the center of the sharp rise of x- Thus we see that 
both the saturation and the increase of x detected by fiSK 
could be qualitatively explained by exchange anisotropy. 

To summarize, our measurements in ZnCu3(OH)gCl2 
reveal an anisotropic intrinsic spin magnetization with 
Xl > xl possibly due to anisotropic g factor. At fields 
above 2 kG a CW temperature can be consistently de- 
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FIG. 6: (Color online) A plot of the magnetization detected 
by muon spin rotation versus temperature (black squares), 
and simulation data for antiferromagnetic kagome lattice with 
Ising-like exchange anisotropy as in Ref. [281 ] (gray lines). In 
the inset the normalized critical temperature versus the ex- 
change anisotropy is shown. 



termined in two different directions. By mean-field ap- 
proximations we were able to show that this phenomenon 
can be explained only by anisotropic super-exchange con- 
stants where J z > Jj_. This anisotropy can explain the 
main features of the susceptibility determined by /iSR. 
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